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Introduction 

Sound  localization  requires  the  computation  of  interaural  time 
differences  (ITDs)  for  frequencies  below  ~1  to  10  kHz.  This  is 
performed  by  binaural  cells  in  the  avian  Nucleus  Laminaris  (NL),  and 
its  mammalian  homologue,  the  Medial  Superior  Olive  (MSO). 

An  “ITD  discriminator”  neuron  should  fire  when  inputs  from  two 
independent  neural  sources  coincide  (or  almost  coincide),  but  not  when 
two  inputs  from  the  same  neural  source  (almost)  coincide.  A  neuron 
that  sums  its  inputs  linearly  would  not  be  able  to  distinguish  between 
these  two  scenarios. 

This  is  a  biophysical  model,  using  the  program  NEURON,  to  examine 
how  NL  neurons  detect  and  report  ITDs,  their  mechanisms  and  their 
limitations. 
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Results 

•  Typical  chick-like  parameters  allow  ITD  discrimination  up  to  2  kHz. 

•  Typical  chick-like  parameters  but  with  barn-owl-like  phase  locking 
allow  ITD  discrimination  up  to  4  kHz  (can  easily  be  pushed  to  6  kHz). 

•  Two  dendritic  non-linearities  aid  ITD  discrimination: 

1)  intra-dendritic  inputs  sum  sub-linearly. 

2)  inter-dendritic  interactions  subtractively  inhibit  out-of-phase  inputs. 

•  Response  to  monaural  input  does  not  require  spontaneous  activity 
from  other  side. 


•  Rate-coded  ITD  tuning  curves  convey  more  information  than  Vector- 
Strength-coded  curves  (despite/due  to  Vector  Strength  enhancement). 
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Model  Description 

The  model  emulates  an  array  of  neurons,  each  with  an  adjustable  number  of  dendrites,  a  soma,  and  an 
axon  with  an  axon  hillock,  a  myelinated  segment,  and  a  node  of  Ranvier.  Each  section  has  an 
adjustable  number  of  equipotential  compartments.  All  geometric,  electrical,  and  channel  parameters 
are  adjustable,  as  are  the  number  of  synapses/dendrite  (~  30),  the  synaptic  locations,  and  the 
distribution  of  synaptic  locations.  Channel  types  include  potassium  (high  and  low  voltage  activated 
[~Kvl.l,  1.2]  and  delayed  rectifier),  sodium,  and  passive.  Values  were  obtained  from  physiological 
studies  of  Nucleus  Magnocellularis  (NM)  and  NL  [Refs.  Rathouz  &  Trussel,  1998,  and  Reyes,  Rubel, 
&  Spain,  1996].  Voltage  dependent  channels  are  specified  by  Hodgkin-Huxley-like  parameters.  Each 
neuron  in  the  array  feeds  into  a  single  inhibitory  neuron,  which  feeds  back  onto  all  neurons  in  the 
array. 

The  stimulus  is  a  pure  tone  of  adjustable  frequency,  with  adjustable  interaural  phase  difference  (or 
contralateral  monaural  stimulus  with  variable  ipsilateral  spontaneous  activity).  More  complex  stimuli 
can  be  easily  introduced. 

The  synapses  fire  with  conductance  proportional  to  an  alpha-function,  with  adjustable  time  constant, 
peak  conductance,  and  reversal  potential.  The  excitatory  synapses  fire  as  individual  Poisson  processes, 
with  probability  rate  given  by  a  modified  sinusoid,  with  adjustable  amplitude  and  vector  strength.  The 
inhibitory  neuron  is  a  simple  integrate- and-fire  type. 


The  implementation  uses  the  program  NEURON  and  has  a  graphical  user  interface  for  controlling  the 
parameters  and  running  the  model.  NEURON  provides  a  real-time  display  of  data  and  analysis 
including  the  potential  at  various  locations,  the  two  stimuli,  the  synaptic  firings,  spike  counters,  period 


histograms  of  synaptic  firings  and  the  action  potentials,  and  their  vector  strengths. 

C/W? 


W  Center  for  Auditory 
and  Acoustic  Research 


Institute  for  Systems  Research 
University  of  Maryland 


NEURON  Panels 


#  dendrites 


in 


■a 


Length  [Den]  (um) 


JF 


'0 


Diameter  [Den]  (um) 


Ax.  Resist.  [Den]  (ohm  cm) 


gL  [Den]  (S/cmA2) 


gK  LVAjn  [Den]  (S/cmA: 


ur 


F 

£pr 


'0 


2UU 


'S 


U.UUUb 


'0 


gK  HVAm  [Den]  (S/cmA2)  Fj]  |U.U3 
Compartments  [Den]  ~|^j  |3U""I 


3d 

ip 


F 


'0 

"a 


lambda  (um)  [ 


'0 


28«.b« 


■0 


#  Ex.  Synapses/dendrite  J  |  |3U 

 0 

Center  [Ex  Syn]  (0,1)  1 _ J  |  O.b 

0 

Distribution  [Ex  Syn]  (0,1)  j| _ |  |1 

0 

tau  [Ex  Syn]  (ms) 

QF 

0 

gmax[ExSyn]  (uS)  |_j  |0.0ib 

 0 

e  [Ex  Syn]  (mV) 

J|-1U 

■0 

Duration  [Ex  Syn]  (ms)  J  | 

0 

Ow? 

W  Center  for  Auditory 
and  Acoustic  Research 


Length  [Hillock]  (um)  I  |  |C 

0 

0 

Diameter  [Hillock]  (um)  Irj 

3 

30 

Ax.  Resist.  [Hillock]  (ohm  cm) 

UpF 

0 

gLeak  [Hillock]  (S/cmA2)  1 

JJ0.0006 

0 

gNa_m  [Hillock]  (S/cmA2)  1 

\\]1H 

30 

gK  m  [Hillock]  (S/cmA2)  j 

|  I0-52 

0 

#  Compartments  [Hillock] 

JF 

30 

Length  [Myelin]  (um)  j  pTJO 

0 

Diameter  [Myelin]  (um)  p  j 

3 

0 

Ax.  Resist.  [Myelin]  (ohm  cm) 

0 

gLeak  [Myelin]  (S/cmA2)  |  ||/.be-0b 

0 

C  [Myelin]  (uF/cmA2)  M  | U.U1 2b 

0 

#  Compartments  [Myelin]  1 _ |  1 1 0 

0 

Length  [Node]  (um)  _ |  |2 

=30 

Diameter  [Node]  (um)  1 

7 

-0 

Ax.  Resist.  [Node]  (ohm  cm) 

p3 

0 

gLeak  [Node]  (S/cmA2)  I 

(U.UUUb 

0 

gNa m  [Node]  (S/cmA2) 

IF55 

0 

gK_m  [Node]  (S/cmA2)  jj 

[TO - 

30 

#  Compartments  [Node]  1 

1133 

0 

Length  [Soma]  (um)  |  |  1 15 

0 

Diameter  [Soma]  (um)  1 

T5 

0 

Ax.  Resist.  [Soma]  (ohm  cm) 

JjjFT 

0 

gK  LVA_m  [Soma]  (S/cmA2) 

Lj|u 

330 

gK  HVA  m  [Soma]  (S/cmA2) 

jj  |0 

330 

gLeak  [Soma]  (S/cmA2)  1 

1 0.0006 

30 

gNa_m  [Soma]  (S/cmA2) 

IF 

30 

gK_m  [Soma]  (S/cmA2)  1 

3 

0 

#  Compartments  [Soma]  1 

13 

0 

eNa  (mV)  1 

TD 

0 

eK  (mV)  [rj 

T5D 

0 

eLeak(mV) 

p6D - 

0 

alphaO  HVA  (/ms 

1  j  |u.n 

0 

alphaVHalf  HVA  (mV)  1 

-19 

— 0 

alphaK  HVA  (mV) 

pF 

0 

betaO  HVA  (/ms) 

J  (U.1U3 

■0 

betaVHalf  HVA  (mV)  M  p 

iy 

30 

betaK  HVA  (mV) 

P  F 

0 

alphaO  LVA  (/ms) 

pP 

'0 

alphaVHalf  LVA  (mV)  1 

^50 - 

0 

alphaK  LVA  (mV) 

P|2i.8 

'0 

betaO  LVA  (/ms) 

J  13 

0 

betaVHalf  LVA  (mV)  M  p 

BU 

3 0 

betaK  LVA  (mV) 

Pi14 

0 

qlOHVA  H 

l 

..0 

TO  HVA  (C)  Fj] 

3 

0 

q1 0  LVA  |[J 

0 

rnFfcrg 

3 

0 

Stimulus  Frequency  (Hz)  jQj  1 1 000 

0 

Stimulus  Phase  Ipsi  (Deg)  II 

|0 

1  ▲  1 

Stimulus  Phase  Contra  (Deg)  1 

JF 

10 

Stimulus  Vector  Strength  (0,1) 

Ul° 

•1  0 

Probability  Rate  (per  ms)  IF]  |0.bb 

=0 

Generic  Parameter 

J  I1 

0 

Action  Pot.  Threshold  (mV)  ||  J  |-us 

0 

Period  Histogram  bins  |_J  1 1 6 

00 

Ignore  spikes  before  (ms) 

|ib 

0 

Cells  per  Array  |_J  |2 

0 

Integration  factor  (uS)  I  1  [ 

UP 


-0 


tau  [In  Syn]  (ms) 


gmax  [In  Syn]  (uS) 
e  [In  Syn]  (mV)  jg[ 


fOU8 


'0 

“0 


'0 


#  Primary  Values  | _ |  1 1 1 

0 

Parameter  to  vary: 

stimFreq 


First  Value  1 _ |  jbbb.bb 

-0 

Last  Value  i  1  |2UUU 

0 

3|  Vary  with  Log  Scale 

Slave  parameter  to  vary: 

--none- 


Slave  First  Value 

 0 

Slave  Last  Value 

JP - 

“ 0 

j]  Vary  with  Log  Scale 

Use  Secondary  in  Parallel 

#  secondary  values  j _ j  1 1 

0 

Parameter  to  vary: 

-none-- 

First  Value  ;_j  |b00 

0 

Last  Value  1  1 12000 

30 

J|  Vary  with  Log  Scale 

Slave  parameter  to  vary: 

-none- 


Slave  First  Value  Ij  1 1 

0 

Slave  Last  Value  ||_|  1 1 

 0 

Vary  with  Log  Scale 

j]  IDen  Follows  stimFreq 

Soma  'g's  from  Axon  &  Dendrite 
a=  VS  Follows  stimFreq  (chick) 

VS  Follows  stimFreq  (owl) 

□  Use  Binaural  Stimulus 


Institute  for  Systems  Research 
University  of  Maryland 


Geometry  &  Connectivity 


A  typical  model  cell  has  2-8  dendrites,  each  20  -  400 
pm  long  and  2-4  pm  in  diameter,  a  spherical  soma  of 
diameter  15  pm,  and  an  axon.  Each  dendrite  has  ~30 
excitatory  synapses.  The  axon  has  an  axon  hillock,  a 
segment  with  myelination,  and  a  node  of  Ranvier.  The 
output  feeds  into  an  integrate- and-fire  inhibitory  cell 
which  feeds  back  to  all  cells  in  the  array. 


Axon  hillock  initiating  spike  from  ipsilateral  current  surge,  despite  contralateral  current  drain 
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out-of-phase  in-phase 


Time  Plots 

A  pair  of  cells  receives  the  same  stimulus  probability  distributions  (here,/=l  kHz). 
The  top  receives  its  inputs  binaurally  in-phase,  and  the  bottom  out-of-phase. 


contra  g  synapse 
ipsi  g  synapse 
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Red  tracks  the  intracellular  potential  in  mid- soma, 
mager  at  the  axon  tip,  and  brown  in  one  dendrite, 
plots  inhibitory  conductance. 

Directly  beneath  are  the  pair  of  presynaptic  stimulus 
probability  distributions.  See  figure  to  right  for  other 
examples  of  stimulus  probability  distributions.  - ► 


The  bottom  8  curves  of  each  graph  show  realized  synaptic 
currents  (note  spread  from  Poisson  distribution). 
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Stimulus  probability  distributions 
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out-of-phase  in-phase 


__  Period  Stimulus  Histograms 

The  same  pair  of  cells  (and  stimulus),  tracked  for  250  ms: 
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Vector  Strength  (VS)  measures  phase  locking,  between  0  and  100%. 

The  in-phase- stimulus  cell  increased  VS  over  its  inputs  (a  little  too  well). 
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Parameter  Space 


The  dimensionality  of  the  parameter  space  is  high,  with  only  a  small  region 
biologically  relevant. 

To  limit  the  search  through  parameter  space,  most  trials  compare  a  pair  of  identical 
cells  receiving  identical  stimuli,  but  one  in-phase  and  the  other  out-of-phase.  Trials 
using  monaural  stimulus  do  not  need  pairs. 


Most  trials  co-vary  the  dendritic  length  and  stimulus  vector  strength  with  the 
stimulus  frequency,  using  experimentally  derived  relations. 


A  linear  fit  between 
lo g(dendritic  length )  and 
log  {best  frequency) 

Data  from  Smith  &  Rubel,  1979. 


A  linear  fit  between 
vector  strength  and 
log  {best  frequency) 


Data  from  Koppl,  1997, 
Warchol  &  Dallos,  1990. 
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ITD  Discrimination 
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Good  ITD  discrimination 
(high  ratio)  until  ~2  kHz. 

At  right,  keeping  all 
parameters  the  same,  but 
using  Barn  Owl  vector 
strength  gives  good  ITD 
discrimination  until  ~4  kHz. 
(One  can  go  up  to  ~6  kHz  250 
simply  by  adding  more 
dendrites.) 
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Non-Linearities 
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Opposite  dendrite’s  effect  is  subtractive 
“ subtracts  when  nothing  positive  to  add ” 
Works  at  all  frequencies,  including  high 
New  result 


Both  effects  prevent  many  inputs  from 
right  side  wrongly  causing  cell  to  fire 
without  inputs  from  the  left  side. 

Reduction  in  “ false  positives  ” 


Synaptic  inputs  add  sub-linearly 
more  inputs  don ’t  add  as  much  you  ’d  think” 
Works  only  at  low-middle  frequencies 
Found  by  Agmon-Snir  et  al. 
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__  Intra-Dendritic  Sub-Linearity  _ 

EPSCs  from  the  excitatory  synapses  sum  This  non-linearity  arises  from  a  low  synaptic 
sub-linearly  on  entering  the  dendrite.  reversal  potential  (EExSyn). 


synapses  i  dendritic 
a  b  |  result 

1  ^  0  ->  1 


0  t  i  ->  1 


1  t  1  ->  1.5 


See  Agmon-Snir,  et.  al.  1998 
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ExSyn 


“(V-E^  Syn) 


Concurrent  EPSCs  raise  the  potential  V, 
diminishing  the  effect  of  each  EPSC.  At  the 
extreme,  if  V  >  EExSyn,  PSCs  are  effectively 
inhibitory. 

The  effect  is  strongest  for  longer  dendrites 
(i.e.  more  electrically  isolated). 


works  best  at  lower  frequencies 


(Synaptic  depression  is  another  example  of  a  sub-linearity.) 
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Sub-Linearity  Results 
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#  Shifting  the  synaptic  reversal  potential  upwards  reduces  the  sub-linearity, 
worsening  the  ratio  of  in-phase/out-of-phase  firing  rates. 

•  The  effect  is  present  only  at  lower  frequencies. 
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Inter-Dendritic  Subtraction 


Current  flowing  into  the  soma  from 
dendritic  EPSPs  is  subtracted  by  the 
opposite  dendrite  (if  without  its  own 
EPSPs).  The  subtraction  is  greater  when 
the  opposite  dendrite  had  recent  EPSPs. 


dendrites  i  somatic 
a  b  |  result 
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0  4  1  ->  0.5 


1  4  1  ->  2 
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This  non-linearity  arises  from  voltage-dependent  K+ 
channels.  When  there  are  no  ESPSs  in  the  opposite 
dendrites,  the  channels  are  somewhat  activated, 
acting  as  a  mild  current  sink.  When  there  were 
recent  ESPSs  in  the  opposite  dendrite,  the  channels 
are  strongly  activated,  acting  as  a  large  current  sink. 

This  gives  a  hierarchy  of  firing  rates 
for  different  stimuli: 
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The  out-of-phase 
rate  is  suppressed 
relative  to  the 
monaural  rate. 


The  opposite 
dendrite  acts  as  a 
current  sink. 


Subtraction  Results 
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Dendrite  removed 

Normal 

High  K+ 

Conductance 


Note  that  the  cells  fire  well  with  no  stimulus  on  the  opposite  side.  The  cell  is  not  just  a 
coincidence  detector,  it  is  an  ITD  discriminator:  it  does  not  need  spontaneous  activity  on  the 
opposite  side  in  order  to  fire. 


The  effect  is  present  at  all  frequencies. 

(The  meeting  of  in-  and  out-  rates  at  ~  2kHz  is  a  consequence  of  poorly  phase-locked  inputs.) 
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Phase  Locking 


The  responses  to  a 
360°  range  of  phase 
difference  give  ITD 
tuning  curves 
(Interaural  Phase 
Difference  + 
frequency  gives  ITD) 

Rate  coding  and 
vector- strength 
coding  give  two 
separate  ITD  tuning 
curves. 

The  “rate”  ITD 
curves  are  more 
sharply  tuned  than 
the  “vector- strength” 
ITD  curves  (note  that 
the  vector  strength  is 
not  reliable  when  the 
firing  rate  is  low). 
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300  - 


Vector  Strength 
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841  Hz  87  um  51%  VS 
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Parameters  used  were  not  tuned  to  give  the  most  accurate  possible  results.  These  results  show 
1)  an  over-enhancement  of  output  VS  over  the  input  VS,  making  the  VS-coded  ITD  tuning 
curves  appear  extra  flat,  and  2)  an  over-suppression  of  rates  for  nearly  out-of-phase  inputs, 
which  makes  the  rate-coded  ITD  tuning  curves  look  extra  sharp  (compared  to  experiment) 
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Vector  Strength  [%] 


Synaptic  Location 


There  are  many 
“experiments”  one  can 
do  with  a  model  that  are 
either  difficult  or 
fundamentally 
impossible  in  the 
laboratory. 

In  this  case  we  compare 
In-Phase/Out-of-Phase 
ratios  for  3  interesting 
but  “impossible” 
scenarios. 

In  all  but  the  black 
(normal)  case,  every 
excitatory  synapse  on 
each  dendrite  has  been 
put  at  the  same  point. 

In  the  3  colored  cases, 
they  are  all  at  the 
dendrite  base,  dendrite 
center,  and  dendrite 
tip. 
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dendritic  length  [fim] 


200  100  50 


Base 

Middle 

Tip 

Normal 


The  only  difference  is  for  long  dendritic  lengths/low  frequencies, 
as  expected,  since  the  intra-dendritic  non-linearity  requires  an 
electrically  isolated  dendrite. 

Base: 

The  intra-dendritic  non-linearity  is  diminished,  the  out-of-phase 
rate  goes  up  (the  number  of  false  positives  goes  up),  and  the  In- 
Phase/Out-of-Phase  ratio  goes  down. 

Middle  &  Tip: 

The  more  isolated  the  synapses  from  the  soma,  the  higher  the 
effect  of  the  intra-dendritic  non-linearity,  and  the  more  the  In/Out 
ratio  goes  up. 
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Dendritic  Length 


The  Intra-dendritic 
sublinearity  leads 
to  an  optimal 
dendritic  length,  as 
shown  by  Agmon- 
Snir  et  al. 

For  every  stimulus 
frequency  there  is 
a  dendritic  length, 
longer  than  which, 
performance  no 
longer  increases. 
The  effect  is  most 
pronounced  at 
lower  frequencies. 
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ITD  Discrimination — Barn  Owl 


dendritic  length  [pm] 


1 00  50  25 


In-Phase  Rate 
Out-Phase  Rate 

log(In/Out) 
In-Phase  VS 


dendritic  length  =  20  |tm,  16  dendrites 


The  top  plot  shows  ITD 
discrimination  using 
chick-like  parameters, 
but  merely  increasing  the 
vector- strength  to  that  of 
the  barn  owl. 


The  right  plot  shows  ITD 
discrimination  going  up 
to  ~6  kHz  by  simply 
adding  more  dendrites. 
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